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Abstract 

Dito and Turrubiates recently introduced an interesting model of the 
dissipative quantum mechanics of a damped harmonic oscillator in phase 
space. Its key ingredient is a non-Hermitian deformation of the Moyal star 
product with the damping constant as deformation parameter. We com- 
pare the Dito- Turrubiates scheme with phase-space quantum mechanics 
(or deformation quantization) based on other star products, and extend 
it to incorporate Wigner functions. The deformed (or damped) star prod- 
uct is related to a complex Hamiltonian, and so necessitates a modified 
equation of motion involving complex conjugation. We find that with this 
change the Wigner function satisfies the classical equation of motion. This 
seems appropriate since non-dissipative systems with quadratic Hamilto- 
nians share this property. 



1 Introduction 



The quantum mechanics of dissipative systems has been studied intensively; for 
reviews, see [H [21 13]. Work started soon after the birth of quantum mechanics 
and continues today. 

Perhaps the most fundamental approach is to consider the system of interest 
as interacting with an appropriate reservoir. Then the well-known quantization 
methods that are valid for non-dissipative, closed systems can be applied to the 
system plus reservoir as a whole. Effective equations of motion for the system 
can then be found by integrating out the reservoir degrees of freedom, while 
making appropriate physical assumptions. 

Another technique is to work backwards, and derive the effective equations 
by adapting quantization procedures to non-dissipative systems. The adapted 
quantization procedures should, in the end, agree with more fundamental treat- 
ments. Provided they do, they would be helpful, as shortcuts for the more 
fundamental derivations, and possibly more. 

Canonical (operator) quantization has been adapted to dissipative systems 
either by using time-dependent Hamiltonians, complex Hamiltonians, or by 
modifying the canonical Poisson brackets, and the corresponding operator com- 
mutatorsjj 

We concern ourselves here, however, with a different quantization method: 
quantum mechanics in phase space, or deformation quantization (see Sect. 2). 
Specifically, we study Dito and Turrubiates [5| recent adaptation of deforma- 
tion quantization to the paradigmatic dissipative system, the damped harmonic 
oscillator. 

The important innovation they introduce is a non-Hermitian 7-deformation 
of the Moyal star product of non-dissipative deformation quantization, where 
7 denotes the damping constant. Their damped star product is built in turn 
on a 7-deformation of the classical Poisson bracket, and so recalls the modified 
brackets of canonical operator quantum mechanics adapted to dissipation that 
were just mentioned. Perhaps this is not surprising, since canonical classical me- 
chanics is a key ingredient of deformation quantization; the algebra of quantum 
observables is described as an ?i-deformation of the classical Poisson algebra on 
phase space. Unlike in other schemes, the one proposed by Dito and Turrubiates 
only modifies the classical bracket and thereby the corresponding quantum star 
product - it uses the undamped Hamiltonian, for example. The Dito-Turribiates 

^ Some relatively recent works with the latter approach are [H [5l l6l Ff). 
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proposal is at least economical, since in other formulations involving modified 
brackets, extra structure must be input. 

This paper is organized as follows. The next section is a brief review of 
ordinary (undamped) quantum mechanics in phase space. Section 3 studies the 
Dito-Turrubiates scheme and extends it to include Wigner functions and their 
evolution. The final section is a discussion, consisting of a concise summary, 
and a list of some of the many questions that remain. 



2 Phase-space quantum mechanics 

We now give a brief review of phase-space quantum mechanics, to establish 
our notation and provide the essential background for next section's discussion 
of the Dito-Turrubiates scheme. Throughout this paper our treatment will be 
limited to one non-relativistic particle travelling on the real line with coordinate 
q, and conjugate momentum p. 

For pedagogical introductions to phase-space quantum mechanics (or defor- 
mation quantization), see [21 [TU], for examples; and for more advanced reviews, 
see the book [11], and references therein, including the pioneering work [12j . 

In our context, operators will be expressible entirely in terms of the operator 
position q and momentum p. These more general operators will also be indicated 
by carets. Phase-space quantum mechanics has a direct relationship with the 
more commonly used operator quantum mechanics. The Wigner transform W 
maps operators to functions and distributions on phase space: 

W(f) = f{p,q) . (1) 



Here / is an operator, and /(p, q) is the corresponding phase-space distribution, 
sometimes called the symbol of /. The most important property of the Wigner 
transform is encoded in 

W(/g) = (W/) * (W.g) . (2) 

Consequently, operators can be Wigner-mapped to symbols whose algebra is 
homomorphic to the original operator algebra, as long as the symbols are mul- 
tiplied using the so-called star-product (^-product) 

* = exp [y(ag9p - 9p5g)] . (3) 

Operators can therefore be replaced by functions/distributions on phase space, 
and quantum mechanics can be done in phase space, without reference to oper- 
ators. 
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The Moyal star product * inherits associativity from the operator product, 
by (HI). Since {fg)^ = P , we also have 



{f * g) g * f ■ (4) 

Here / indicates the complex conjugate of /. We say that * is Hermitian. 

The inverse of the Wigner transform is familiar in operator quantization. 
The so-called Weyl map W^^, satisfies 

>V-i(/ * 5) - • (5) 

The image of a polynomial in q and p of this Weyl map is the corresponding 
Weyl-ordered operator. The famous Moyal star product * of ([3]) and ([5]) is 
therefore intimately related to Weyl ordering. With other operator orderings, 
other ^-products arise. 

The operator encoding information about the quantum state is the density 
matrix, and its Wigner transform is the central object of phase-space quantum 
mechanics. Up to normalization, it is the celebrated Wigner function. 

To consider evolution, start with the Liouville Theorem for the classical 
phase-space density pc'- 

_ = _ + = 0. (6) 

The Dirac correspondence then yields the equation of motion of the density 
operator (matrix) p 

+ [p,H] = 0. (7) 
The Wigner transform of this is the evolution equation for the Wigner function 

p = p{p, q; t) 

ih^ + [p, i?], - , (8) 
where the ^-commutator is 

[a , b]^ = a * b — b * a . (9) 

Equivalently, the phase-space version of the Dirac quantization rule, 

K.l^Ktk, (10, 

leads directly from to ([8|). 
Defining 

ad^.M.9 if , 9]* , (H) 
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([5]) can be rewritten as 

= ad^HP i^^P- (12) 

The solution is just 

p{p,q;t) = expj -^ad^H I p(p,g;0) = exp { t £ } g; 0) . (13) 

Using the associativity of *-niuhiphcation, this reduces to 

p{p,q;t) = U{p,q;t) * p{p,q;0) * U{p,q;-t) , (14) 

with 

m...,^t^X^y ^e.,m(^). (15, 

Here 

ExpM(a) (16) 

n=0 



denotes the so-called ^-exponential. Note that U{p, q; t) — U{p, q\ —i), and 



U{p,q-t)*U{p,q-t) = U{p,q;t)*Uip,q;t) - 1. (17) 

Alternatively, substituting ([M]) into ([8]) yields the dynamical equation 

ihdtU{p,q;t) = H * U{p,q;t) , (18) 



and its complex conjugate. Its solution, the symbol of the propagator, is (|T5|1 . 
The spectrum of energies and the Wigner functions of stationary states can then 
be found by the spectral decomposition (or Fourier-Dirichlet expansion) 

U{p,q-t) = Y,p^{p,q)e-^^"^ , (19) 

E 

where 

H * p^ = * H = Ep^ . (20) 

For the simple harmonic oscillator, with Hamiltonian 

1 

H = — + -muj^q^ , (21) 

the well-known results are 

X 1 f 2iJtan(wi/2) 1 

U{p,q-t) = exp ' > , 22 

cos(a;t/2) |^ inuj J 
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and 

"^^ ''•'{^) '''^{-'^) ■ <^^' 

where En = hu}{n + 1/2), and L„ is the n-th Laguerre polynomiaL 

The time evohition invokes more than the pure-state, or diagonal Wigner 
functions, however. We need the off-diagonal 

satisfying the *-eigen equations 

H*Pe,e' = EPe.e'^ Pe,e''^H = E'p^^^,; (25) 

so that ^ ~ Pe' Assuming that the *-eigen functions are complete, we 
expand 

p{p.q;t) = ^e.e'^^'^Pe,e'^p^i) ■ (26) 

E.E' 

Substituting into the equation of motion then gives the time-evolved Wigner 
function 



9; = X! ^E E' (^) '^^P 



E,E' 



-i{E~ E')t 



Pe.e'^P^'i)- (27) 



To close this section, we will discuss alternatives to the Moyal star product 
(see [TU] and [H], e.g.). First, consider operator-ordering ambiguities. Since 
the Moyal product is intimately related to Weyl ordering, a different operator 
ordering will lead to a different star product. For example, if we use the so- 
called standard ordering, where every q is placed to the left of every p, we find 
instead the standard star product 

Since the two orderings can be simply related, there is a nice relation between 
the two star products: 

T,{f *,9) - (T^/) * (T^g) , (29) 
involving the invertible transition operator 

= exp|-^a,9p| . (30) 

An ordering change preserves a bi-grading in the powers of p and q. By the 
Heisenberg commutation relation, h has bi-grade (1,1) and so that of is 
(0,0). 



6 



This kind of connection between different star products can be extended. In 
general, a star product S is considered equivalent to the Moyal one *, if we have 

f{f *9) = {ff)~*{fg) . (31) 

Here T stands for an invertible transition operator. If ?, like *, is an h- 
deformation of the point-wise product of functions, then 

f = 1 + 0{n) . (32) 

Associativity of ? is guaranteed by the relation ([3T|) . if T is invertible: 

figih = f((f-V) * (f-'g) * (T-i/i)) . (33) 

If T is real, then S will be Hermitian, as * is. * „ is an example of a non-Hermitian 
product, since is not real. 

By the Liebniz product rule, the equivalence relation pip can be rewritten 

as 

s = * f-^[d] f[d + d] f-^[d] . (34) 

Here d stands for either dq or dp, 

f[d + d] :=f 1 ^ ^ , (35) 

and similarly for r^^[(5] and T^^[d]- For example, 

T^[d + d] = cxp I (dg + dg) (dp + 9p) I , (36) 

by dSOl). 

Eqn. ([M]) makes it clear that equivalent star products * and * only differ by 
a factor that is left-right symmetric. Consequently, it can be shown that 

limMi ^ l^Ml = (37) 
h^o in h^o in 

the Poisson bracket. 

The equivalence between S and * defined by (j3ip is known as c-equi valence. 
It is important because the Moyal ^-product ^ is the unique, associative de- 
formation of the point- wise multiplication of functions on (and R^"), up to 
c-equivalence. 

The c in c-equivalence stands for cohomological. As such, it is a mathemat- 
ical, rather than a physical equivalence. For example, consider the Husimi star 
product *^ (see [13], e.g.) 

*H = *exp||(s2 dqdq +^ dpdp)^ , (38) 
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related to * by the transition operator 



H 



exp 



(39) 



It is also easy to see that *^ satisfies the condition ([57)1 . But *^ is used with 
the Husimi phase space distribution 



(40) 



which does not encode precisely the same physics as the Wigner function p{p, q; t). 
Eqn. ((40|) can be rewritten as 



1 

■kTi 



dp'dq' p{p',q';t) 



X exp 



1 

h 



[q - q'f 



s^{p-p' 



(41) 



indicating that the Husimi distribution is a smoothed version of the Wigner 
function, coarse-grained by a squeezecH Gaussian weighting in phase space. So, 
although *^ is c-cquivalent to the Moyal product, it describes different physics. 



3 The damped harmonic oscillator in phase- 
space quantum mechanics 

We start in this section with a review of the Dito and Turrubiates [S] model 
of damping in a quantum harmonic oscillator in phase space, with comments 
added. A new contribution will then be described: our incorporation of Wigner 
functions and their evolution. First, a naive proposal will be examined, and 
rejected, since it leads to an unacceptable evolution equation. The equation 
of motion will then be modified, and explained. Finally, its consequences are 
outlined. One such consequence is that the Wigner function of the damped 
harmonic oscillator follows the canonical flow. This property is also common to 
all non-dissipative systems having quadratic Hamiltonians, including the simple 
(undamped) harmonic oscillator. Since the damped and undamped harmonic 
oscillator and such systems are treated similarly in some other approaches, this 
seems physically reasonable. 

Dito and Turrubiates [H] describe the damped harmonic oscillator using the 
Hamiltonian of the undamped simple harmonic oscillator. In a sense then, they 

More precisely, the distribution functions for general s were introduced in I14II15I . while 
the Husimi distribution has s = 1. We call the Gaussian weighting of II41I I squeezed because 
it is proportional to the Wigner transform of a squeezed state. 
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describe its dissipation as kinematics. More precisely, the dissipation is encoded 
in a deformation of the classical Poisson bracket of the harmonic oscillator and 
its consequent quantum ^-product. 

The initial observation is that the classical equations of motion of the damped 
harmonic oscillator can be written as 

9 = H}^ = p/m ; p ^ {p, H}^ = -mu'^q - 2jp , (42) 

using the undamped Hamiltonian (|2ip . The price paid is that we must use a 
deformed Poisson bracket 

{fi9}-y ■= {/,.9} - 2-fmdpfdpg 
= f (^dq dp - dp dq - dp dp^ g , (43) 

with the damping constant 7 as deformation parameter. Notice that the de- 
formed bracket is no longer skewH and doesn't obey the Jacobi identity. 

This result applies to a general class of classical systems: if we use any 
Hamiltonian of the form 

H = ^ + Vi,) , (44) 

with an arbitrary potential V{q), the same damping term 2m^q is still the only 
modification of the equation of motion for q{t). We will restrict attention here 
to the simple harmonic oscillator and its damped version, however. 

Just as the Moyal product is related to the Poisson bracket, the Dito- 
Turrubiates damped ^-product is obtained by exponentiation of the deformed 
classical bracket (l43l): 



exp 



Yidqdp - dpdq -27m dpdp)] = H.e-*''^'"^-^- . (45) 



Why exponentiate? In the undamped case, it is necessary for the homomorphism 
of the *-algebra with the operator algebra. Also, the Moyal ^-product is the 
unique associative deformation of the point-wise multiplication of functions on 
R^", up to isomorphism (by transition operators - see below). No similar result 
for the damped product is known; neither is the operator algebra for the damped 
case. Exponentiation is therefore an assumption. 



^ The deformed bracket can be obtained from the Poisson bracket by the substitutions 

dq^dq — 2am7 dp and dq^dq +2l3m-f dp, for any real a and /3 such that Q + /3 = 1. The 
sign difference in front of 7 between the substitutions for the left- and right-acting derivatives 
foreshadows our proposal 17011 . 
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On the other hand, suppose that we posit 

*7 := p(^Y^dqdp - dpdq -2jm dpdp)^ , (46) 

for some function F. Since {•, -j-y {■,■} as 7 — > 0, requiring that hm-y^o *7 = 
* selects the exponential function. The possibihty that the form is too 
restrictive remains, however. 

One key result of [5] is 

T(/ * g) = {Tf)*,{Tg) , (47) 

with 

3 = exp 



( = 5.^). (48) 



This can be proved easily using That is, the Dito-Turrubiates damped star 
product is c-equivalent to the Moyal star product. As discussed in the previous 
section, c-equivalence does not imply physical equivalence, and so the damped 
star product has the potential to describe the dissipation of a damped harmonic 
oscillator. 

So, what effect does the Dito-Turrubiates transition operator T have? First, 
it is clear that T does more than change the ordering, since it has no fixed 
bi-grade in the powers of p and q. Its effect is therefore more profound. An 
important example is 

+ V{q)) = ^ + V{q) - ^^ ; (49) 
\ 2m / 2m 2 

(see also ((55)) )H| T transforms a real Hamiltonian into a complex one. 

For such a conversion to be possible, it is necessary, but not sufficient, that 
T =/= T. Another consequence of a non-real transition operator is that is not 
Hermitian: 



a^jb = b *_-y a ^ b *j a . (50) 

Its c-equivalence with the Moyal star product ensures that the damped star 
product *j is associative, for any value of the damping constant 7. For example, 
will be useful later because of ([50]) . and it is also associative. However, 
problems exist if any two different damping parameters 7 and 7' are used. Care 
must be taken with regard to the order of multiplications. Clearly, 

/*7(1*7'.9) / (/*7l)*7'ff' (51) 



* This last result is also interesting in its own right, since it recalls a different approach to 
damped systems that uses complex Hamiltonians (see [ll[2], and also |16p. 
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if 7 7^ 7', for example. A unique product such as a^^b c does not (automat- 
ically) exist. 

Use of a naive substitution 

KH.-^^ (52) 

leads to 

= ^ + l[p,,i^k, (53) 

where indicates the Wigncr function of the damped harmonic oscillator. This 
equation of motion integrates to 

Pi{p,q\t) = U^{p,q;t)*.yPj{p,q;0)*jU^{p,q;-t) . (54) 

Here 

[/,(p, 9; t) ^ _ = Exp[*,] (55) 

satisfies the damped analogue of i.e. 

z7i9t [/-y {p,q]t) = H U-y (p, g; i) . (56) 

Dito-Turrubiates solve ([55]) , without discussing the evolution of Wigner func- 
tions. They then Fourier-Dirichlet expand the solution to find the spectrum of 
eigenvalues and eigenfunctions of 

Using (|47| . one finds 

r 1 f -itH\ ( , , /~itT-^(H) 
U^{p,q;t) = ExpK]f^— j = T^ExpMf 

= t(ExpM (zM^±iM) ) . eW.T(c/(p,,;0) ■ (58) 

Using (HH) and (gS]), this gives [5] 

exp(7i/2) 



UJp,q;t) 



cos 



(wt/2) [l + ^tan(wt/2)] 

X exp<| l^^^^M^ ( + 2 ) ^ (59) 

^ I 2m[l + ^tan(tjt/2)] 2 ^ ' < ^ ' 



The expansion of Uj gives the solutions of (|57p in terms of the simple harmonic 
oscillator analogues: 

P7.« = t(p^) , (60) 
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with ^ ^ 

E^,n = En + = - [{2n + l)u + Z7] . (61) 

For example, for n = Dito and Turrubiates |:8^ find 
2 

^^•y.o ^1 - 2i-i/uj 

X expl --^ ^ ^ — -- + imt^vl 1 ■ (62) 

^\ ?iw V2rn[l - 2i7/w] 2 ^/J ^ 

However, there is a fundamental problem with the basic equation of motion, 
eqn. (j53p . For the damped harmonic oscillator it yields 

° = + ^ + z^t^^'^l* + *7?^5,9,p^ . (63) 

Notice that the term proportional to the damping constant 7 is not real. But 
the eigenvalues of the density matrix are the populations, and they must be real. 
Time evolution must preserve the reality of the density matrix. The imaginary 
damping term in (|63p means that a real density matrix does not remain real as 
it evolves. 

An equally important flaw is revealed by the limit 

lim = {H,p^} + (64) 

This is a direct consequence of the c-equivalence of and *. But this shows 
that there is no damping in the classical limit. The classical limit of (|55|) is not 
correct! 

To try to understand better the origin of the difficulties with the purported 
equation of motion let us consider how it might be "derived". Notice that 
(|53p would result from the undamped evolution equation ([5]) by the substitution 

^ = T{p) , (65) 

a generalization of (|60p . This simple identification is appealing in part because it 
is similar to (|40p . that defining the Husimi phase-space distribution p^{p,q;t). 
The Husimi equation of motion is found directly from that for the Wigner 
function, by substituting (j40p . Additional terms arise in the Husimi equation 
of motion compared to that for the Wigner function, since the coarse-graining 
evolves in time |13| . Just as the "twisting" by Th modifies the equation of 
motion, so does application of the Dito- Turrubiates transition operator T. The 
physical damping is meant to be introduced that way. 
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Following (|¥T|) . one might hope for an interpretation of T{p) as the Wigner 
distribution coarse-grained in momentum space. That point of view is not 
sensible, however, since the weighting would have to be Gaussian-like with an 
imaginary exponent. This does point to the origin of the main problem, however: 
unlike the Husimi transition operator Th, the Dito-Turrubiates T is not real. 

As a consequence, eqn. (1491) can hold: T transforms a real Hamiltonian into 
a complex one. That is the crucial point here, we believe. Compare to a similar 
situation - if a non-self-adjoint Hamiltonian is used, the equation of motion of 
the density operator is modifiedH to 

= ih^ + pH - p . (66) 

By analogy, we should consider]^ 



= Pi*^H - p^^^H . (67) 
Notice that the right-hand side of this equation is purely imaginary, so that 

^ = -|lni(p,*.i/ ) . (68) 
This implies that the reality of the Wigner function 

Pi ^ (69) 
is preserved in evolution^ We can therefore also write 

dp-y 

- ^ p^*iH - H p (70) 



With this prescription it is easy to show that the classical limit makes sense for 
the simple harmonic oscillator Hamiltonian: 

p *iH-H*^^p {p H}.^ - {H,p }_.^ 

tl Jn = 2 = ^^-^ ■ 

We emphasize that the relation is then not obeyed. 

Let us now attempt to argue for ([70]) from other grounds. As a starting 
point, let us assume that the Liouvillc Theorem still holds, and try to modify 
the argument that led to ^ to justify the damped equation of motion ([70]) . An 



^ See |17) . for example, where the analogous modified equation of motion for a Heisenberg 
operator was shown to lead to a quantum anomaly. 

^ This form may possibly be related to the bi-orthogonal quantum mechanics discussed by 
Curtright and Mezincescu [ 18] . For related work in phase space, see [19] and 1201. 

^ According to ([68}, if p had a non-zero imaginary part, it would not evolve. 
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important advantage of the Dito-Turrubiates method is that it only modifies the 
classical brackets. That advantage would be lost if we need to input something 
to replace the Liouville Theorem. Therefore, we assume the equation of motion 
is 

+ {p^^cH}^ = . (72) 



dt dt 

The phase-space version of the Dirac quantization rule pI7|) above should there- 
fore be deformed to 



{a,b}, ^ - = - (73) 



for real observables a and b, in order to recover the equation of motion ((67l 
One might be tempted to write 



P-y{p,q;t) = U^{p,q;t)*-^ p^{p,q-0) U^{p,q;t) (74) 

as a real solution to the equation of motion (j70p . But this expression is ambigu- 
ous at best, as shown by the discussion around eqn. (ISTj) . 

Luckily, however, a formal solution to (j70p can be written, by deforming the 
undamped solution of (|llj|13p . If we define 

ad^^'[/]ff / <? - 5*7/ , (75) 

then ^ is 



= ad^j\H]p^ =: thC^ p^ . (76) 



The solution is just 



p^{p,q;t) = exp|^adJ-^)[H]| p^(p,9;0) = exp {t C^} p^{p,q;0) . (77) 

There is no associative ambiguity in this explicit solution. 

A more explicit result can be given immediately for the damped harmonic 
oscillator, having 



2 9 P f d r, „ 



The quantum evolution is ([70]) . but that reduces to ([7^ . with p^ ,, p^. The 
quantum Wigner function satisfies the classical equation of motion. The Wigner 
function at time t is therefore simply 

P^{p,q;t) = f{Pc{-t),qc{-t)) , (79) 
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where Pc{t) and qdt) characterize the classical (daraped) trajectories in phase 
space. Thus the quantum damped harmonic oscillator follows the classical back- 
ward flow of the phase-space coordinates. Perhaps this is not surprising, since 
for any quadratic Hamiltonian, such as that of the undamped simple harmonic 
oscillator, the same results holds (see [TT], e.g.). 

4 Discussion 

We first summarize. We have shown how to describe the dynamics of Wigner 
functions in the Dito-Turrubiates scheme [5]. The non-Hermitian damped star 
product *j is c-equivalent to the Moyal product *. Therefore, if the evolution 
equation of the damped Wigner function only involves ^-^-commutators, only 
Poisson brackets survive in the classical limit, rather than the damped bracket 
{•, - j-y. The classical limit would then be damping-free, and therefore incorrect. 

However, the damped transition operator T transforms the simple harmonic 
oscillator Hamiltonian into a complex one, according to Consequently, the 
evolution equation does not involve ^^-commutation, but must instead be (j70p . 
This ensures that a real Wigner function remains real as it evolves, and the 
classical limit is correct. Not only is the classical limit correct, it is exact. The 
damped harmonic oscillator in this formalism therefore follows the classical flow, 
a property shared with non-dissipative systems in phase space having quadratic 
Hamiltonians. 

Most helpful to us were (i) comparisons of the damped star product *^ 
with other physical star products that are also c-equivalent to the Moyal *, viz. 
the standard product and the Husimi product and (ii) the Heisenberg 
equation of motion for an operator observable modified for a non-Hermitian 
Hamiltonian (see [T7], e.g.). 

Let us now conclude with a few of the many questions that remain. We hope 
progress can be made toward their answers. 

Should the damped *-y-product be obtained by exponentiating ih{-, ■}j/2? In 
the case at hand, much of the structure of the 7-deformed star product is irrel- 
evant. Since we use the quadratic Hamiltonian (PT|) . the only terms that enter 
are those that are up to quadratic in d and d- In this sense, our main result has 
a certain robustness. More work on this question would be helpful, however. 

Is there a deformed structure analogous to the Heisenberg- Weyl group that 
is relevant to the damped case? In the undamped case, the Heisenberg- Weyl 
group is the structure of paramount importance. The Moyal ^-product simply 
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provides a ^-realization of that group. As discussed around ([?T|) . associativity 
can be a problem in the damped case. But if (j77p is an appropriate guide, 
perhaps 



exp 



\^8id[^'>[ap + bq] } 



^cp-\-dq 



^cp-\-dq —ih(ad—bc-\'2m'yac) 



(80) 



can serve the purpose. Here a, b, c and d are constants, so that e"P+'"' and e'^^^'^'^ 
♦-represent elements of the Heisenberg-Weyl group. When 7 — » in a form 
of the defining ^-relations of the Heisenberg-Weyl group is recovered. 

Can a solution to the damped equation of motion {7C^ analogous to the un- 
damped formula \2'T^ be written ? It is clear from (j77p that the *-eigen equation 
of ad^^''''' [iJ] is relevant, so one can start there. One possible ansatz follows, 
although it may only have relevance for very small 7. Suppose we find the 
time-independent off-diagonal Wigner matrix elements p 
satisfying 



i;£,£' 



H p 



!■£,£' 



^ f^r,£,£' ' 



1\£:£ 



,*jH = S' p , 



f;£,£' 



where the complex eigenvalue^ are 

£ = E + iX , £' = E' + i\' 



(81) 



(82) 



with E, E' G M, and A, A' G M+. If the real eigenfunctions of ([8T|l exist and are 
complete at all times t, then we can expand 



(83) 



£,S' 



The equation of motion then yields 



£,£' 



A{£-£')t 



5]i?^_^,(0)exp 



-i(E-E')t 





■-(A + A')i" 


exp 





(84) 



If this conjecture is correct, then the physically important *-eigen equations are 
those given in ([5T|) . It would not be clear then that the phase-space functions 
of (|60p are directly relevant. 

The ansatz of (1811184^ reduces to the undamped system when h ^ Q. It is 
also devoid of unphysical stationary states - the imaginary parts of eigenvalues A 



* For the role such complex eigenvalues play in another description of the damped harmonic 
oscillator, see |21| . Wigner functions involving such eigenvalues are considered in I22| . 
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and A' both contribute to the damping of the dynamics. Furthermore, the system 
is consistent with ([7^ if {£ — £')/fL is independent of h. Such eigenvalues were 
obtained in as indicated in (pT|) . However, the corresponding solutions (|60p 
are not real - see (|^. e.g. 

Can the Dito-Turrubiates scheme he related to other quantization methods 
for dissipative systems (see m\^\^)? Certainly, hints of connections with other 
models are apparent. For example, Dekker's use of a complex Hamiltonian is 
recalled by T{H) = H — ihj/2. Bateman's doubled system of a damped and 
anti-damped oscillator comes to mind from a,dSj\H]p ~ H*^p — p *^-y H; in 
this last expression, however, the anti-damped (7 — > —7) system is dual rather 
than extra/ auxiliary. The importance of resonances has been emphasized by 
[22I [2T] and others, and seems relevant to the *~^- and *_.^-eigen equations of 
the previous paragraph. Finally, if the correspondence (j65p were correct, then 
we would be able to define a damped Wigner transform 

:= T-^W = Wf-\ (85) 

where 

= exp|-^(ad[x])'}. (86) 

The form of this T appears to be consistent with Tarasov's 23 proposal to 
include superoperators in the operator formulation of quantum mechanics in 
order to adapt it to dissipative systems. 

Can other dissipative physical systems be treated in a similar way? One very 
useful example might be spin systems, with relaxation. 

Instead of dissipation, might other physical effects, such as decoherence, be 
describable in the Dito-Turrubiates manner? 
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